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In condensed matter physics gauge symmetries other than the U(l) of electromagnetism are of 
an emergent nature. Two emergence mechanisms for gauge symmetry are well established: the 
way these arise in Kramers- Wannier type local-global dualities, and as a way to encode local 
constraints encountered in (doped) Mott insulators. We demonstrate that these gauge structures 
are closely related, and appear as counterparts in either the canonical or field-theoretical language. 
The restoration of symmetry in a disorder phase transition is due to having the original local 
variables subjected to a coherent superposition of all possible topological defect configurations, with 
the effect that correlation functions are no longer well-defined. This is completely equivalent to 
assigning gauge freedom to those variables. Two cases are considered explicitly: the well-known 
vortex duality in bosonic Mott insulators serves to illustrate the principle. The acquired wisdoms 
are then applied to the less familiar context of dualities in quantum elasticity, where we elucidate the 
relation between the quantum nematic and linearized gravity. We reflect on some deeper implications 
for the emergence of gauge symmetry in general. 
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I. INTRODUCTION. 

Breaking symmetry is easy but making symmetry is 
hard: this wisdom applies to global symmetry but not 
to local symmetry. The study of systems controlled by 
emergent gauge symmetry has become a mainstream in 
modern condensed matter physics. Although one dis- 
cerns only the gauge invariance of electromagnetism in 
the ultraviolet of condensed matter physics, it is now very 
well understood that in a variety of circumstances gauge 
symmetries that do not exist on the microscopic scale 
control the highly collective physics on the macroscopic 
scale. An intriguing but unresolved issue is whether the 
gauge structures involved in the Standard Model of high 
energy physics and perhaps even general relativity could 
be of such an emergent kind. 

The mechanisms for gauge symmetry emergence fall 
into two broad categories: (1) the "stay-at-home" gauge 
invariance associated with (doped) Mott insulators and 
the gauge fields associated with the slave-particle theories 
encountered in this context, describing the fractionaliza- 
tion of the quantum numbers of the microscopic degrees 
of freedom, and (2) the global-to-local symmetry cor- 
respondence encountered in the strong-weak (Kramers- 
Wannier, S-) dualities as of relevance to the quantum 
field theories describing the collective quantum physics 
of condensed matter systems. In the common perception 
these appear as quite different. The purpose of this note 
is to clarify that at least in the context of bosonic physics 
they are actually closely related. In fact, these high- 
light complementary aspects of the vacuum structure, 
and it is just pending whether one views the vacuum ei- 
ther using the canonical/Hamiltonian language (stay-at- 
home) or field-theoretical/Lagrangian (local-global du- 



ality) language. The case is very simple and we will 
illustrate it in section [TTJ with the most primitive of all 
many-particle systems governed by continuous symme- 
try: the Bose-Hubbard model at zero chemical poten- 
tial, or alternatively the Abelian-Higgs duality associ- 
ated with complex-scalar field theory. Although we do 
not claim any new result in this particular context, the 
freedom to switch back and forth between the Lagrangian 
and Hamiltonian viewpoints yields some entertaining vis- 
tas on this well-understood theory. To make the case that 
it can yield new insight, we apply it in section [TIT] to the 
less familiar context of dualities in quantum elasticity. 
This deals with the description of quantum liquid crys- 
tals in terms of dual condensates formed from the trans- 
lational topological defects (dislocations) associated with 
the fully ordered crystal. Using the Lagrangian language 
it was argued that such quantum nematics are equivalent 
to (linearized) Einstein gravitj^. Here we will demon- 
strate that this is indeed controlled by the local symme- 
try associated with linearized gravity: translations are 
gauged, turning into infinitesimal Einstein transforma- 
tions. 

This could have been a very short communication 
but we wish to address a readership with diverse back- 
grounds. We ther efore first review at length the bas ics of 
Mottness (section II A ) and vortex duality (section II B ) 



which should be quite familiar for the condensed mat- 



ter physicists, presenting in section II C our argument 
revealing how the stay-at-home gauge is encoded in the 
vortex duality. Section [III| is devoted to duality in quan- 
tum elasticity and the relation to gravity. We again take 
the time to review the basics since this subject is un- 
familiar even for the condensed matter mainstream. In 
section |II C| we review the basics of these duality struc- 
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tures while in section |IIIB| we turn to the way this is 
related to gravitational physics, employing the insights 
derived from vortex duality to show why the effective 
space realized in the quantum nematic behaves like the 
spacetime of general relativity when it is nearly flat. 

II. VORTEX DUALITY VERSUS BOSE-MOTT 
INSULATORS. 

In this section we will first review two standard views 
on the physics of the Bose-Hubbard modeP, with the in- 
tention to confuse the reader. We first highlight the emer- 
gence of the compact U(l) stay-at-home gauge which 
emerges in the Bose-Mott insulator in a trivial way when 
one sticks to the canonical, second quantized language 
focusing on the large-?/ limit. In the second subsection 
we turn to the field-theoretical Abelian-Higgs duality, re- 
viewing the standard argument that in 2+1D the quan- 
tum disordered partner of the superfluid is actually dual 
to a superconducting condensate formed from the vor- 
tices of the superfluid, which interact via effective C/(l) 
gauge fields- 8 . These should be equivalent descriptions 
of the Bose-Mott insulator, but at face value this is far 
from obvious. In the last subsection we present a very 
simple resolution of this conundrum, having some inter- 
esting ramifications for the way one should think in gen- 
eral about "Mott insulators" in field theory. 

A. The stay-at-home gauge of the Bose-Mott 
insulator. 

A mainstream of the gauge theories in condensed mat- 
ter physics dates back to the late 1980s when the commu- 
nity was struggling with the fundamentals of the problem 
of high-T c superconductivity. It was recognized early on 
that this has to do with doping the parent Mott insu- 
lators and this revived the interest in the physics of the 
Mott insulating state itsefP^. The point of departure 
is the Hubbard model for electrons, 

H F n = -t (CAv + c] a c la ) + Uj2h tt n li , (1) 

<ij>a i 

describing fermions cj CT on site i with spin a, hop- 
ping on a lattice with rate t, subjected to a strong local 
Coulomb interaction U. Here hi a = cj^Cio- is the fermion 
number operator. When U 3> t and when there is on 
average one fermion per site, it is easy to recognize the 
Mott insulator: a fermion will localize at every site, while 
it costs an energy ~ U to move a charge which will there- 
fore not happen at macroscopic scales. This simplicity is 
deceptive: viewed from a general perspective the effect 
of this "projective renormalization" is profound. Due to 
the dynamics as dominated by the strong local repul- 
sions the particle number becomes locally conserved. In 
the large- U limit, 



x;^i*)=i*>. (2) 

a 

where | ^) is the ground state wave function. Although 
the on-site fermion number is not a conserved quantity of 
the microscopic Hamiltonian Eq. ([I]) it becomes locally 
conserved at low energy. This truly emergent local con- 
straint/conservation law can be imposed by a compact 
[7(1) gauge field a ia , 

c f ct e iai " 

hi = ^ ^L^icr h t . (3) 

a 

In this effortless way one discovers that a gauge sym- 
metry emerges that controls the physics at long distances, 
while it is non existent at the microscopic scale. This is 
the point of departure of a mainstream school of thought 
in condensed matter physics. There is still a dynamical 
spin system at work at low energies. Using various "slave- 
constructions" it was subsequently argued that quantum 
spin liquids characterized by fractionalized excitations 
can be realized when the resulting compact U(l) gauge 
theory would end up in a deconfining regime. 

A much simpler problem is the Bose-Hubbard model 
describing spinless bosons created by 6] hopping on a 
lattice with a rate t subjected to an on-site repulsion U, 

HBH = -tJ2 h lbj+Uj2^- ( 4 ) 

<ij> i 

Again fii — b\bi is the boson number operator. We as- 
sume in the remainder that the system is at "zero chem- 
ical potential", meaning that on average there is an in- 
teger number of bosons no per site. Recently this prob- 
lem attracted attention in the context of cold bosonic 
atoms on a optical lattice^, while it describes equally 
well Josephson junction networks^]. In the large-?/ limit 
one obtains the same simple picture for the Bose-Mott 
insulator as for its fermionic sibling, except that the spin- 
less Bose variety is seen as completely featureless since it 
does not seem to break a manifest symmetry while low 
energy degrees of freedom are absent. However, it does 
share the trait with the fermionic Mott insulator that 
when U/t 1 a low energy compact U(l) gauge invari- 
ance is generated since the number of bosons per site is 
sharply quantized, 

b\ -> b\e ia * , 

h % = b\b t ->• hi . (5) 

In 2+1 and higher dimensions it is well established 
that the model Eq. Q is characterized by a quantum 
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phase transition as function of the coupling g = U /not 
at a critical coupling g c . The Mott insulator is estab- 
lished for g > g c while for small couplings g < g c the 
zero temperature state is a superfluid. In this regime it 
is convenient to rewrite the Hamiltonian in phase rep- 
resentation, b\ = |6|e 4( ^, to obtain the phase dynamics 
model (J — not), 



Ha, — J 



<i]> 



COS 



(6) 



subjected to the quantization condition [(pi, hi] = iSij. 
For g < g c the global £7(1) symmetry associated with the 
phase 4> will break spontaneously and this results in the 
superfluid. In terms of the original boson the symmetry 
breaking implies that it develops a vacuum expectation 
value (VEV), 



n e 



i<t>o 



(7) 



where 4>o the phase shared by the condensate, and no 
the average density. Comparing with Eq. ([5]), it is ob- 
vious that the global U(l) symmetry that breaks spon- 
taneously in the phase condensate acquires an emergent 
compact U(l) gauge status in the "number condensate" 
realized in the Mott insulator. This Bose-Hubbard model 
reveals perhaps the most primitive mechanism for the 
emergence of a local symmetry from microscopic global 
symmetry. 



B. Phase dynamics duality: the superconductor as 
the dual of the superfluid. 

Let us now turn to the held-theoretical formulation of 
the same problem. The quantum partition sum is given 
in terms of the Euclidean path integral, 



Z 



V<t> e 



i drd d x c 



(8) 



The Lagrangian C is derived from Eq. ^ departing 
from the superfluid phase by realizing that in the Legen- 
dre transformation hi — > id T (j)i, followed by a naive coarse 
graining. After scaling out the phase velocity c p h ~ y/UJ 
and defining the coupling constant as g = U/ J, one finds, 



C = kg ' 



(9) 



in a rclativistic short hand notation (x, fj, — t,x, y, ...), 
while the field <f> is compact with periodicity 2tt. This 
theory just describes the Goldstone boson/phase mode 
of the ordered superfluid. Recently we showed that this 
theory is subjected to a global-local duality in all di- 
mensions equal to and larger than 2+irJ^. Here we will 



just consider the familiar 2+1D "Abelian-Higgs" or "vor- 
tex" duality case since it can be regarded in the present 
context as fully representative for the higher dimensional 
cases as well. One first transforms Eqs. ( 8|9 | to vortex 
coordinates. After a Hubbard-Stratonovich transforma- 
tion, C = + iJfid^cf). The auxiliary fields J M just 
represent the supercurrents dual to the phase field. This 
is easily seen by factorizing the phase field in smooth 
and multivalued configurations, <f> — 4>sm + 4>m\- The 
smooth (j) sm enters as a Lagrange multiplier which can be 
integrated out yielding the conservation law d„Ju, = 0, 
which is just the supercurrent continuity equation. In 
2+1D this can be imposed by parametrizing the cur- 
rent in terms of a non-compact U(l) gauge field as 
J a = e fj,vxdi/A\. Inserting this in the action, one finds 
after some straightforward reshufflings, 



r = i. F f + iA T v 



d v d. 



AfMV 



(10) 



where F M „ = d^Ay — d^A^, the field strength of A^ 
which is now sourced by vortex currents; one recognizes 
that the non-integrability in the phase field just cor- 
responds with the vorticity. 

We encounter here yet another emergent gauge theory, 
but this one is completely unrelated to the gauging of 
the phase field in the canonical formulation. We have 
just written down an action describing how the topo- 
logical excitations of the superfluid (vortices) mutually 
interact. Since the phase field deformation due to the 
vortex is long ranged, the vortex-vortex interaction is 
long ranged as well. The above derivation shows that 
these long range vortex interactions are actually indistin- 
guishable from the electromagnetic interactions between 
electrically charged particles in 2+ ID. It is in a way just 
a convenience that these can be parametrized in terms 
of gauge potentials, and these non-compact U(l) fields 
have clearly no relation to the emergent compact U(l) 
fields imposing number quantization in the Mott insulat- 
ing phase. 

The next step in vortex duality is to consider what 
happens when the coupling constant g increases. The 
vortices are the unique sources of quantum fluctuations 
and for small couplings they occur as small vortex-anti- 
vortex loops in Euclidean spacetime. For growing cou- 
pling constant these loops will become larger until at g c a 
"loop blowout" occurs where they become infinitely large. 
This is the quantum phase transition, where at larger 
couplings a "tangle of free (anti-)vortex world lines" is 
found in spacetime. This Bose condensate of vortices 
just corresponds with the quantum disordered phase, and 
because of the gauge interactions between the vortices 
this is just the same entity as a condensate of rclativistic 
bosons interacting via electromagnetic gauge fields. This 
condensate is governed by the Ginzburg-Landau- Wilson 
action, 
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Sv = I dr I dx 



\{d„ + iA^ Y \ 2 + m^l^vl 2 + ^v|*v| 4 + jiF^F, 



r 



(ii) 



where Vty is the vortex condensate field, while the 
gauge fields descend from Eq. (10). This is in turn 
just a relativistic superconductor or "Higgs phase" . In 
the condensate the vortex-matter field can be integrated 
out and the effective action for the massive photons can 
be written in a gauge invariant fashion in terms of the 
original supercurrents J M as, 



Mott 



dr I dx 



(12) 



Describing a doublet of propagating photons with a 
mass mMott- The conclusion is that in 2+1D a quan- 
tum disordered superfluid (in terms of phase) can be 
equally well be interpreted as a relativistic superconduc- 
tor formed from the vortices. 



C. The vortex condensate as the generator of the 
hidden stay-at-home gauge. 



with the vortices which are in turn in a perfect non-local 
relation with the original phase variables. However, we 
know precisely what this dual superconductor is in terms 
of those phase variables. We can resort to a first quan- 
tized, world line description of the vortex superconduc- 
tor, putting back "by hand" the phase variables. This 
constitutes a tangle of world lines of vortices, warping 
the original phases, and eventually we can even map that 
back to a first quantized wave function written as a co- 
herent superposition of configurations of the phase field. 
To accomplish this in full one needs big computers^!, 
but for the purposes of scale and symmetry analysis the 
outcomes are obvious. 

The penetration depth Ay of the dual vortex super- 
conductor is just coincident with the typical distance be- 
tween vortices. At distances much shorter than Ay the 
vortices do not scramble the relations between the phases 
at spatially separated points and at these scales the sys- 
tem behaves as the ordered superfluid, 



Up to this point we have just collected and reviewed 
some well-known results on phase dynamics. However, 
at first sight it might appear as if the matters discussed 
in the two previous subsections are completely unrelated. 
Departing from the Bose-Hubbard model the considera- 
tions of the previous subsection leave no doubt that in 
one or the other way the dual vortex superconductor can 
be adiabatically continued all the way to the strongly 
coupled Bose-Mott insulator of the first subsection. The 
standard way to argue this is by referral to the excitation 
spectrum. The Bose-Mott insulator is characterized by 
a mass gap ~ U (at strong coupling), and a doublet of 
"holon" (vacancy) and "doublon" (doubly-occupied site) 
propagating excitations being degenerate at zero chemi- 
cal potential. The vortex superconductor is a relativistic 
U{l)/U(l) Higgs condensate in 2+1 dimensions charac- 
terized by a Higgs mass (a gap) above which one finds a 
doublet of spin- 1 (left and right helical) "vector bosons". 
In this regard there is a precise match. However, in the 
canonical formalism one also discovers the emergent U(l) 
invariance associated with the sharp quantization of lo- 
cal number density in the Mott insulator. What has hap- 
pened with this important symmetry principle in the vor- 
tex superconductor? 

The answer is: the emergent compact U(l) gauge sym- 
metry of the Mott insulator is actually a generic part of 
the physics of the relativistic superconductor. 

The argument is exceedingly simple. The stay-at-home 
gauge does not show up explicitly in the Higgsed action 
describing the dual vortex condensate, for the simple rea- 
son that all the quantities in this action are associated 



(b\r)b(0)) -t constant , r < \\ 



(13) 



However, at distances of order Ay and larger, the vac- 
uum turns into a coherent quantum superposition of 
"Schrodinger cat states" where there is either none, or 
one, or whatever number of vortices in between the points 
and r where we ask the question of how to correlate the 
phases of bosons, see Fig. [T] We have arrived at expos- 
ing the simple principle which is the central result of this 
paper: since the vortex configurations are in coherent su- 
perposition, the phases acquire a full compact U(l) gauge 
invariance. The computation is already done using the 
canonical formalism, and all that remains is to under- 
stand the physical concept: focus on the direction of the 
phase at the origin and look at the phase arrow at some 
distance point r. Consider a particular configuration of 
the vortices, and in this realization the distant phase will 
point in some definite direction which will be different 
from the phase at the origin as determined by the par- 
ticular vortex configuration. However, since all different 
vortex configurations are in coherent superposition and 
therefore "equally true at the same time" , all orientations 
of the phase at point r are also "equally true at the same 
time" and this is just the precise way to formulate that 
a compact U(l) gauge symmetry associated with (j> has 
emerged at distances Ay. 

The implication is that via Eq. ([5| the emerging stay- 
at-home gauge invariance implies that in the Higgs con- 
densate the number density associated with the bosons 
condensing in the dual superfluid becomes locally con- 
served on the scale Ay. The Mottness therefore sets in 
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FIG. 1: In the vortex condensate the correlation of the phase between a point A and another point B a distance r apart is in 
a superposition of have zero, one or any number of vortices in between. As such the phase at B with respect to that at A is 
completely undefined: it has acquired a complete gauge invariance in the sense that any addition to phase is an equally valid 
answer 



only at scales larger than this Ay. Notice that this mech- 
anism does in fact not need a lattice: it is just a generic 
property of the field theory itself, which is independent 
of regularization. In fact, the seemingly all important 
role of the lattice in the standard reasoning in condensed 
matter when dealing with these issues is a bit of tunnel 
vision. It focuses on the strong-coupling limit where for 
large U, Ay — > a, the lattice constant. However, upon 
decreasing the coupling the stay-at-home gauge emerges 
at an increasingly longer length scale Ay, to eventually 
diverge at the quantum phase transition. Close to the 
quantum critical point the theory has essentially forgot- 
ten about the presence of the lattice, just remembering 
that it wants to conserve number locally which is the 
general criterium to call something an insulator. In fact, 
Mottness can exist without a lattice altogether. A rel- 
ativistic superconductor living in a perfect 2+1D con- 
tinuum is physically reasonable. Since duality works in 
both directions, this can be in turn viewed as a quantum 
disordered superfluid, where the number density associ- 
ated with the bosons comprising the superfluid becomes 
locally conserved. 

By inspecting closely this simple vortex duality we 
have discovered a principle which might be formulated 
in full generality as: the coherent superposition of the 
disorder operators associated with the condensation of the 
disorder fields has the automatic consequence that the or- 
der fields acquire a gauge invariance associated with the 
local quantization of the operators conjugate to the oper- 
ators condensing in the order field theory. We suspect 
that this principle might be of use also in the context of 
dualities involving more complex field theories. To sub- 
stantiate this claim, let us now inspect a more involved 
duality which is encountered in quantum elasticity, where 
the principle reveals the precise reasons for why quantum 
liquid crystals have dealings with general relativity. 



III. QUANTUM NEMATIC CRYSTALS AND 
THE EMERGENCE OF LINEARIZED GRAVITY. 

Einstein himself already forwarded the metaphor that 
the spacetime of general relativity is like an elastic 
medium. Is there a more literal truth behind this 
metaphor? In recent years the mathematical physi- 
cist Hagen Kleinert has been forwarding the view that 
quite deep analogies exist between plastic media (solids 



with topological defects) and Einsteinian spacetim d 16 * 17 [ 
There appears room for the possibility that at the Planck 
scale an exotic "solid" (the "world crystal") is present, 
turning after coarse graining into the spacetime that we 
experience. 

It turns out that this subject matter has some bearing 
on a much more practical question: what is the gen- 
eral nature of the quantum hydrodynamics and rigidity 
of quantum liquid crystals? Quantum liquid crystals^ 
are just the zero temperature versions of the classical 
liquid crystals found in computer displays. These are 
substances characterized by a partial breaking of spatial 
symmetries, while the zero temperature versions are at 
the same time quantum liquids. Very recently indications 
have been found for variety of such quantum liquid crys- 
tals in experimentfiSHSl] i n the present context we are es- 
pecially interested in the "quantum smectics" and "quan- 
tum nematics" found in high-T c cuprates^ * 22 * 25 * which 
appear to be also superconductors at zero temperature. 
Such matter should be at least in the long-wavelength 
limit governed by a bosonic field theory, and this "the- 
ory of quantum elasticity"" 2 ^"^ is characterized by du- 
alities that are richer, but eventually closely related to 
the duality discussed in the previous section. Departing 
from the quantum crystal, the topological agents which 
are responsible for the restoration of symmetry are the 
dislocations and disclinations. The disclinations restore 
the rotational symmetry and the topological criterium for 
liquid crystalline order is that these continue to be mas- 
sive excitations. The dislocations restore translational 
symmetry, and these are in crucial regards similar to the 
vortices of the previous section. In direct analogy with 
the Mott insulator being a vortex superconductor, the 
superconducting smectics and nematics can be univer- 
sally viewed as dual "stress superconductors" associated 
with Bose condensates of quantum dislocations. 

Using the geometrical correspondences of Kleiner^ 16 * 17 !, 
arguments were forwarded suggesting that the Lorentz- 
invariant version of the superfluid nematic in 2+ ID is 
characterized by a low energy dynamics that is the same 
as at least linearized gravity^. Very recently it was 
pointed o ut tha t this appears also to be the case in the 
3+1D cas <j 30 * 31 l A caveat is that Lorentz invariance is 
badly broken in the liquid crystals as realized in con- 
densed matter physics. This changes the rules drasti- 
cally and although the consequences are well understood 
in 2+lD^3Ml it remains to be clarified what this means 
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for the 3+1D condensed matter quantum liquid crystals. 
The unresolved issue is how the gravitons of the 3+1D 
relativistic case imprint on the collective modes of the 
non-relativistic systems. 

Here we want to focus on perhaps the most fundamen- 
tal question one can ask in this context: although gen- 
eral relativity is not a Yang-Mills theory, it is uniquely 
associated with the gauge symmetry of general covari- 
ance or diffeomorphism. Quite generally, attempts to 
identify "analogue" or "emergent" gravity in condensed 
matter systems have been haunted by the problem that 
general covariance is quite unnatural in this context. The 
gravity analogues in so far identified in condensed matter 
like systems usually get so far that a non-trivial geomet- 
rical parallel transport of the matter is identified, that 
occurs in a "fixed frame" or "preferred metric'''**"™, or 
either this issue of the mechanism of emerging general 
covariance is just not addressecPsHUl, As we will discuss, 
crystals are manifestly non-diffeomorphic. However, the 
relativistic quantum nematics appear to be dynamically 
similar to Einsteinian spacetime. For this to be true, it 
has to be that in one or the other way general covariance 
is emerging in such systems. How does this work? 

In close parallel with the vortex duality "toy model" 
of the previous section, we will explicitly demonstrate in 
this section that indeed general covariance is dynamically 
generated as an emergent IR symmetry. However, there 
is a glass ceiling: the geometry is only partially gauged. 
Only the infinitesimal "Einstein" translations fall prey 
to an emergent gauge invariance while the Lorentz trans- 
formations remain in a fixed frame. This prohibits the 
identification of black holes and so forth, but this symme- 
try structure turns out to be coincident with the gauge 
fix that is underlying linearized gravity. The conclusion 
is that relativistic quantum nematics form a medium that 
supports gravitons, but nothing else than gravitons. 

For this demonstration we have to rely on the detour 
for the identification of the local symmetry generation 
as introduced in the previous section. Different from the 
Bose-Mott insulator, there is no formulation available for 
the quantum nematic in terms of a simple Hamiltonian 
where one can directly read off the equivalent of the stay- 
at-home gauge symmetry. We have therefore to find the 
origin of the gauging of the Einstein translations in the 
physics of the dislocation Bose condensate, but this will 
turn out to be a remarkably simple and elegant affair. 

The remainder of this section is organized as the pre- 
vious one. In section ITlIAl we will first collect the various 
bits and pieces: a sketch of the way that "dislocation du- 
ality" associates the relativistic quantum nematic state 
with a crystal that is destroyed by a Bose condensate of 
dislocations. In section [ill B| we will subsequently review 
Kleinert's "dictionary" relating quantum elasticity and 
Einsteinian geometry, while at the end of this subsection 
we present the mechanism of gauging Einstein transla- 
tions by the dislocation condensate. For simplicity we 
will focus on the 2+1D case; the generalities we address 
here apply equally well to the richer 3+1D case. 





(a)Edge dislocation 



(b)Screw dislocation 



FIG. 2: Dislocation lines (red spheres) in the relativistic 3D 
"world crystal" (two space and one time direction), formed 
by insertion of a half-plane of particles. Shown in red is the 
contour that measures the mismatch quantized in the Burgers 
vector (red arrow) . If the Burgers vector is orthogonal to the 
dislocation line it is an edge location; if the Burgers vector 
is parallel it is a screw dislocation. In non-relativistic 2+ ID 
there are only edge dislocations, since the Burgers vector is 
always purely spatial. 



A. Duality in quantum elasticity: the quantum 
nematic as a dislocation Bose condensate. 

^26H29l 



Let us first introduce the field-theoretical sic 
The theory of quantum elasticity is just the 19th cen- 
tury theory of elasticity but now embedded in the Eu- 
clidean spacetime of thermal quantum field theory. To 
keep matters as simple as possible we limit ourselves to 
the Lorentz-invariant "world crystal" , just amounting to 
the statement that we are dealing with a 2+ ID elastic 
medium being isotropic, both in space and time direc- 
tions, 



Z = 

S e l = 

where, 



Dw e 2 S <= 



dr / dx 



, (14) 



w 



(15) 



are the strain fields associated with the displacements 
u v of the "world crystal atoms" relative to their equi- 
librium positions. As before g is the coupling constant, 
while /i and A are the shear modulus and the Lame con- 
stant of the world crystal, respectively. At first view 
this looks like a straightforward tensorial generalization 
of the scalar field theory of the previous section. For the 
construction of the nematics one can indeed think about 
the displacements as "scalar fields with flavors" since this 
only involves the "Abelian sector" of the theory associ- 
ated with translations. One should keep in mind however 
that one is breaking Euclidean space down to a lattice 
subgroup and this is associated with non-Abelian, infi- 
nite and semi-direct symmetry structure: the full theory 
beyond the dislocation duality is a much more compli- 
cated affair. 
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FIG. 3: 90° disclination in a square lattice. A wedge is in- 
serted into a cut in the lattice. There is now one lattice point 
with five instead of four neighboring sites (red); going in con- 
tour around this point will result in an additional 90° rotation. 
The associated topological charge is the Frank vector, orthog- 
onal to the plane and of size 90°. As the dislocation, in 2+ ID 
spacetime the disclination point will trace out a world line. 



one direction a superfluid and in the other still a solid. To 
complete this outline, when the coupling constant is fur- 
ther increased there is yet another quantum phase tran- 
sition associated with the proliferation of disclinations 
turning the system into an isotropic superfluid. 

Let us now review the "dislocation duality": in close 
analogy with vortex duality, this shows how crystals and 
liquid crystals are related via a weak-strong duality. The 
requirement that disclinations have to be kept out of the 
vacuum is actually a greatly simplifying factor. One 
follows the same dualization procedure for the disloca- 
tions as for the vortices. Hence, we introduce Hubbard- 
Stratonovich auxiliary tensor fields a^ v , rewriting the ac- 
tion as, 



These complications become manifest considering the 
topological defects: the dislocations and disclinations. 
The dislocation is the topological defect associated with 
the restoration of the translations. The edge dislocation 
can be viewed as the insertion of a half-plane of extra 
atoms terminating at the dislocation core. One imme- 
diately infers that it carries a vectorial topogical charge: 
the Burgers vector indexed according to the Miller in- 
dices of the crystal. In 2+1 dimensions the dislocation is 
a particle (like the vortex) and as an extra complication 
the Burgers vector can either lie perpendicular ("edge 
dislocation", Fig. 2(a)) or parallel ("screw dislocation", 
Fig. 2(b) I to the propagation direction of its world line. 
The disclination is on the other hand associated with the 
restoration of the rotational symmetry. This can be ob- 
tained by the Volterra construction: cut the solid, insert 
a wedge and glue together the sides (see Fig. |3|. This 
carries also a vectorial charge (the Frank vector) . Finally 
dislocations and disclinations are not independent. On 
the one hand, the disclination can be viewed as a stack of 
dislocations with parallel Burgers vectors, while the dis- 
location can be viewed as a disclination-antidisclination 
pair displaced by a lattice constant. 

Dislocations and disclinations have however a different 
identity and this makes possible to give a tight, topologi- 
cal definitions of quantum smectic and nematic order. A 
state where dislocations have spontaneously proliferated 
and condensed, while the disclinations are still massive, 
is a quantum liquid crystal. Since a disclination is coinci- 
dent with a "uniform magnetization" of Burgers vectors, 
one cannot have a net density of parallel Burgers vec- 
tors as long as disclinations are suppressed. The Burgers 
vectors of the dislocations in the condensate have to be 
anti-parallel and therefore the dislocation breaks orienta- 
tions rather than rotations, with the ramification that the 
order parameter is a director instead of a vector. Finally, 
when all orientations of the Burgers vectors are popu- 
lated equally in the condensate one deals with a nematic 
breaking only space rotations. When only a particular 
Burgers vector orientation is populated one is dealing 
with a smectic because the translations are only restored 
in the direction of the Burgers vector: the system is in 



S = I drdx 



1 

4^ 



2 



(16) 

where v = A/2(A + fi) is the Poisson ratio. We divide 
the displacement fields (having the same status as the 
phase field in vortex duality) in smooth and multivalued 



parts Up 



.smooth 



,MV 



, and integrating out the smooth 



strains yields a Bianchi identity, 



= 



(17) 



The physical meaning of is that they are the stress 
fields, which are conserved in the absence of external 
stresses as in Eq. (17): the above is just the stress- 
strain duality of elasticity theory. One now wants to 
parametrize the stress fields in terms of a gauge field. 
Since the stress tensor is symmetric this is most nat- 
urally accomplished in terms of Kleinert's double curl 
gauge fields, 



7 in/ f/ifcA^i/K'A 



>d K d K ,B 



(18) 



while the B's are symmetric tensors, otherwise trans- 
forming as 1-form ?7(1) gauge fields. 

To maintain the analogy with the vortex duality as 
tightly as possible, one can as well parametrize it in a nor- 
mal gauge field, a^ v = e^^xd^A 1 ^ with the requirement 
that one has to impose the symmetry of the stress tensor 
explicitly by Lagrange multipliers. Using this route one 
finds that the multivalued strains turn into a source term 
iA v u jy where, 

fj, fj, 



MV 



(19) 



This is just like a vortex current carrying an extra "fla- 
vor" v. This is the dislocation current, where the fla- 
vor indicates the d+1 components of the Burgers vector. 
As the vortices, dislocations have long-range interactions 
which are parametrized by the gauge fields A (or B), 
with the special effect that these are only active in the 
directions of the Burgers vectors. 



The double curl gauge fields have the advantage that 
the symmetry is automatically built in while the "ex- 
tra derivatives" make possible to identify the disclination 
currents. One finds, 



the Franck vector component. The fact that the disclina- 
tion current has "one derivative less" than the dislocation 
current actually implies that disclinations are in the solid 
confined — in the solid, a disclination is like a quark. 



S 



drdx 



4//. 



(20) 

where the "stress gauge fields" B are sourced by a total 
"defect current", 



a o MV 



(21) 



where # M „ is the disclination current, and v refers to 



One now associates a much larger core energy to the 
disclinations than to the dislocations, and upon increas- 
ing the coupling constant a loop blowout transition will 
occur involving only the dislocation world lines — it is ob- 
vious from the single curl gauge field formulation that dis- 
locations are just like vortices carrying an extra "Burgers 
flavor" . To obtain the quantum nematic one populates 
all Burgers directions equally and after some straight- 
forward algebra one obtains the effective action for the 
"Higgsed stress photons" having the same status as Eq. 
(12 1 for the Mott insulator, 



S 



drdx 



'02' 



1 

4^ 



(22) 



where a should be expressed in the double curl gauge field 
B^ v according to Eq. ( 18 ). In terms of the regular gauge 



fields A v , the first term represents a Higgs mass, while 
the second term is like a Maxwell term. Nevertheless, in 
the nematic the disclinations still act as sources coupling 
to the double curl gauge fields. 

Ignoring the disclinations, one finds in 2+1 D that Eq. 



(22 1 describes a state is quite similar to a Mott insulator: 



all excitations are massive, and one finds now a triplet of 
massive "photons". These are counted as follows: there 
are two propagating (longitudinal and transversal acous- 
tic) phonons of the background world crystal, turning 
into "stress photons" after dualization and acquiring a 
mass in the nematic. In addition, the dislocation con- 
densate adds one longitudinal stress photon. 

As it turns out, the rules change drastically upon 
breaking the Lorentz invariance. In a crystal formed 
from material bosons, displacements in the time direc- 
tion u T are absent, and this has among others the conse- 
quence that the dislocation condensate does not couple to 
compressional stress. Instead of the incompressible na- 
ture of the relativistic state, one finds now two massles 
modes in the quantum nematic: a rotational Goldstone 
boson associated with the restoration of the broken rota- 
tional symmetry, and a massless sound mode which can 
be shown to be just the zero sound mode of the super- 
fluid. The non-relativistic quantum liquid crystals are 
automatically superfluids as well and their relation to 
gravity is obscured. 

Turning to the 3+ ID case one finds as extra compli- 
cation that dislocations turn into strings and one has to 
address the fact that the "stress superconductor" is now 
associated with a condensate of strings. One meets the 



same complication in the case of the vortex duality and 
we showed recently how to handle this"". The outcome is 
actually quite straightforward: the effective London ac- 
tions of the type Eqs. ( 12|22 ) have the same form regard- 
less whether one deals with particle or string condensates, 
and these enter through the Higgs term ~ a 2 /d 2 . 

How to interpret the 2+1D relativistic quantum ne- 
matic? There are no low energy excitations and it only 
reacts to disclinations. It has actually precisely the same 
status as a flat Einsteinian spacetime in 2+1D that only 
feels the infinitesimal vibrations associated with gravita- 
tional events far away. Similarly, using the general rel- 
ativity (GR) technology of th e nex t section, it is also 
straightforward to demonstrat e 1 30 * 31 ^ that in 3+ ID one 
ends up with two massless spin-2 modes: the gravitons. 
To prove that it is precisely linearized gravity, let us con- 
sider next the rules of Kleinert that allow to explicitly 
relate these matters to gravitational physics. 



B. The field theory of quantum elasticity and 
geometry: the Kleinert rules. 

Elaborating on a old tradition in "mathematical metal- 
lurgy" , Kleinert identified an intriguing portfolio of gen- 
eral correspondences between the field theory describing 
elastic media and the geometrical notions underlying gen- 
eral relativity. In order to appreciate what comes, we 
need to familiarize the reader with some of the entries of 
this dictionary. For an exhaustive exposi tion w e refer the 
reader to Kleinert's books on the subjectf^ 1 ^. 

GR is a geometrical theory which departs from a met- 
ric gp.ii, such that an infinitisimal distance is measured 
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through, 



ds — Q i±ii dx ^ dx 1/ , 



(23) 



One now insists that the physics is invariant under lo- 
cal coordinate transformations (general covariance) x^ — ¥ 
£,ii{xv)] infintesimal transformations then are like gauge 
transformations of the metric, 



9^v -t 9^ + + d v ^ = g^ + h 



(24) 



Only quantities are allowed in the theory which are in- 
variant under these transformations and insisting on the 
minimal number of gradients, one is led to the Einstein- 
Hilbcrt action governing spacetime, 



S = — 



2k 



d d xdt Rx 



(25) 



where g = det g^ and R the Ricci scalar, while n is set 
by Newton's constant. Together with the part describing 
the matter fields, the Einstein equations follow from the 
saddle points of this action. 

How to relate this to solids? Imagine that one lives 
inside a solid and all one can do to measure distances is 
to keep track how one jumps from unit cell to unit cell. In 
this way one can define a metric "internal" to the solid, 
and the interesting question becomes: what is the fate 
of the diffeomorphisms ("diffs") Eq. |pi|? In order to 



change the metric one has to displace the atoms and this 
means that one has to strain the crystal, 



(26) 



But the strain fields are surely not gauge fields: the 
elastic energy Eq. ( 14 ) explicitly depends on the strain. 



Obviously, the crystal is non-diffeomorphic and it is char- 
acterized by a "preferred" or "fixed" frame. This is the 
deep reason that normal crystals have nothing to do with 
GR. 

In standard GR the objects that are invariant keep 
track of curvature and these appear in the form of curva- 
ture tensors in the Einstein equations. Linearizing these, 
assuming only infinitisimal diffs as in Eq. (24), one finds 



for the Einstein tensor appearing in the Einstein equa- 
tions, say in the 2+1D case to avoid superfluous labels, 



G 



(27) 



One compares this with the disclination current Eq. 



(21) and one discovers that these are the same expres- 



sions after associating the strains u> M „ with the infinitesi- 
mal diffs h^u- This is actually no wonder: at stake is that 
the property of curvature is independent of the gauge 
choice for the metric. One can visualize the curved man- 
ifold in a particular gauge fix, and this is equivalent to 



the fixed frame. The issue is that curvature continues to 
exist when one lets loose the metric in the gauge volume. 

What is the meaning of the dislocation tensor? Cartan 
pointed out to Einstein that his theory was geometrically 
incomplete: one has to allow also for the property of 
torsion. It turns out that torsion is "Cartan-Einstein" 
GR sourced by spin currents and the effects of it turn 
out to be too weak to be observed (see e.g. Ref. H3"|) . 
In the present context, the torsion tensor appearing in 
the equations of motions precisely corresponds with the 
dislocation currents. With regard to these topological 
aspects, crystals and spacetime are remarkably similar. 

However, given the lack of general covariance the dy- 
namical properties of spacetime and crystals are entirely 
different. For obvious reasons, spacetime does not know 
about phonons while crystals do not know about gravi- 
tons, let alone about black holes. A way to understand 
why things go so wrong is to realize that the disclinations 
encode for curvature, while gravitons can be viewed as 
infinitisemal curvature fluctuations. As we already ex- 
plained, disclinations are confined in crystals meaning 
that it costs infinite energy to create curvature fluctua- 
tions in normal solids. 

Let us now turn to the relativistic quantum nematics: 
here the situation looks much better. Gravity in 2+1D 
is incompressible in the sense that the constraints do not 
permit massless propagating modes, the gravitons. We 
also found out that disclinations are now deconfined and 
they appear as sources in the effective action Eq. (22): 
this substance knows about curvature. In fact, one can 
apply similar considerations to the 3+1D case, where two 
massless spin-2 modes are present. The relativistic quan- 
tum nematic in 3+1D behaves quite like spacetime! 

To make the identification even more precise, one no- 
tices that the expression for the linearized Einstein ten- 
sor Eq. (27) is coincident with the expression for the 



stress tensor in terms of the double curl gauge field B^ Vl 
Eq. (18). But now one is dealing with gauge invariance 



both of B^ v and h^v while they are both symmetric ten- 
sors. At least on the linearized level the stress tensor 
is the Einstein tensor. It is now easy to show that the 
Higgs term in the theory of the nematic when expressed 
in terms of the linearized Einstein tensors, 



-> R 



(28) 



actually reduces to the Ricci tensor R, demonstrating 
that one recovers the Einstein-Hilbert action at distances 
large compared to the Higgs scale. Again, this only holds 
in the linearized theory. This works in the same way in 
3+1 (and higher) dimensions which is the ea sy wa y to 
demonstrate that gravitons have to be presentP212il. At 
least the linearized version of the Einstein-Hilbert action 
appears to be precisely coincident with the effective field 
theory describing the collective behavior of the quantum 
nematic! 
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Although this all looks convincing there is still a gap 
in the conceptual understanding of what has happened 
with the geometry of the crystal in the presence of con- 
densed dislocations. The emergence of gravity requires 
that the original spacetime defined by the crystal has to 
become diffcomorphic. The fields as of relevance to the 
dynamics of the nematic are healthy in this regard but 
they belong to the dual side. The analogy with the Mott 
insulator is now helpful: to demonstrate that gravity has 
emerged requires the demonstration that the spacetime 
of the original crystal is diffeomorphic and that is equiv- 
alent to demonstrating that in the vortex condensate the 
superfluid phase acquires a compact U(l) gauge invari- 
ance. The diffeomorphic nature of the stress gauge fields 
telling about the excitations of the quantum nematic has 
in turn the same status as the gauge fields that render 
the vortex condensate to be a superconductor. 

The good news is that we can use the same "first quan- 
tization" trick that helped us to understand the emer- 
gence of the stay-at-home gauge in the vortex condensate 
to close this conceptual gap. As for the vortices, it is 
easy to picture what happens to the metric of the crystal 
when the coherent superpositions of dislocation config- 
urations associated with the dual stress superconductor 
are present. Let us repeat the exercise at the end of 
the previous section, by comparing how two points some 
distance apart communicate with each other, but now fo- 
cusing on the metric properties. This is illustrated in Fig. 
|4j imagine that no dislocation is present between the two 
points and one needs N jumps to get from one point to 
the other. However, this configuration is at energies less 
than the Higgs mass of the quantum nematic necessar- 
ily in coherent superposition with a configuration where 
a dislocation has moved through the line connecting the 
two points: one now needs N + 1 hops and since these 
configurations are in coherent superposition "N = N+V 1 
and the geometry is now truly diffeomorphic! 

However, there is one last caveat. Although transla- 
tional symmetry is restored in the quantum nematic, the 
rotations are still in a fixed frame and even spontaneously 
broken! This is different from full Einstein gravity: in 
real spacetime also the Lorentz transformations (rota- 
tions in our Euclidean setting) are fully gauged. In order 
to understand this point, let us start from special rel- 
ativity, which has the global symmetry of the Poincare 
group comprising translations and Lorentz transforma- 
tions. The translations form a subgroup, such that trans- 
lational and rotational symmetry are easily distinguish- 
able. More precisely, the generators of translations are 
ordinary derivatives which commute [8^,8^] = 0. In 
many ways, going from special to general relativity is 
from going from global to local Poincare symmetry^. 
Indeed, referring to elasticity language, it seems to make 
sense to restore first translational and then rotational 
symmetry, ending up in a perfectly locally symmetric 
"liquid" state. 

However, it has long been known that such "gauging of 
spacetime symmetry" is very intricate, which has to do 



with the definition of locality under such transformations. 
What happens is that local coordinate tranformations of 
the form — > ^(x v ), which are in fact local transla- 
tions, can also correspond to local rotations. The local 
translations no longer form a subgroup, as the generators 
of translations should be augmented to those of parallel 
translations, defined bjS^, 

D (l = d lt +T l f x f KX , (29) 

where is the connection and f K \ is the generator of 
local rotations. Such modified derivatives do not com- 
mute, and two consecutive translations may result in a 
finite rotation. Such symmetry structure is actually at 
the heart at everything non-linear happening in Einstein 
theory including black holes. 

Going back to what we now know of the quantum ne- 
matic, it is clear that it cannot correspond to full GR, 
since rotational symmetry as reflected by disclinations 
is still gapped. Nevertheless, the identification between 
quantum nematics and linearized gravity is in perfect 
shape. Linearized gravity is a special and somewhat 
pathological limit of full GR, as it only applies to nearly 
globally Lorentz symmetric systems. It was quite some 
time ago realized that such systems are symmetric under 
global Lorentz transformations and infinitesimal coordi- 
nate transformations (see ch. 18,35 in Ref. PfS"|) . This is 
equivalent to fixing the Lorentz frame globally yet allow 
for infinitesimal Einstein translations. Under such condi- 
tions the equations of motion of linearized gravity follow 
automatically. 

Here we have demonstrated that linearized gravity — 
a very peculiar limiting case of GR — is actually literally 
realized in a quantum nematic. The deeper reason is 
that in a quantum nematic the rotational symmetry of 
(Euclidean) spacetime is global and even spontaneously 
broken, while the restoration of the translational sym- 
metry by the dislocation condensate has caused the fixed 
frame internal coordinate system of the crystal to turn 
into a geometry that is characterized by a covariance ex- 
clusively associated with infinitesimal translational coor- 
dinate transformations. 



IV. CONCLUSIONS. 

In so far as the Abelian-Higgs (or "vortex" ) duality is 
concerned we have presented here no more than a clar- 
ification. Living on the "dual side", where the Bosc- 
Mott insulator appears as just a relativistic supercon- 
ductor formed from vortices, the emerging stay-at-home 
local charge conservation from the canonical representa- 
tion in terms of the Mott insulating phase of the Bose- 
Hubbard model is not manifestly recognizable. However, 
the dual vortex language contains all the information re- 
quired to reconstruct precisely the nature of the field con- 
figurations of the "original" superfluid phase fields which 
are realized in the vortex superconductor. By inspecting 
these we identified a very simple but intriguing principle. 
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FIG. 4: In the dislocation condensate (quantum nematic), the distance between two points (green dots) is in a coherent 
superposition of having zero, one or any number of half-line insertions (light blue) or dislocations (red dot) in between them, 
and therefore the number of lattice spacings in between them is undefined. This is equivalent to having the Einstein translations 
fully gauged: there is a diffeomorphism between configurations with any number of lattice spacings in between the two points. 



The local charge conservation of the Mott insulator, as- 
sociated with the emergent stay-at-home compact £7(1) 
gauge symmetry, is generated in the vortex condensate by 
the quantum mechanical principle that states in coherent 
superposition "are equally true at the same time" — the 
Schrodinger cat motive. 

We find this simple insight useful since it yields a some- 
what more general view on the nature of strong-weak du- 
alities. We already emphasized that Mott insulators as 
defined through the local conservation of charge do not 
necessarily need a lattice. One does not have to dig deep 
to find an example: our dual superconductor is just a rel- 
ativistic superconductor in 2+ ID, which is in turn dual 
to a Coulomb phase that can also be seen as a supcrfluid. 
The charge associated with this supcrfluid is locally con- 
served in the superconductor, regardless of whether the 
superconductor lives on a lattice or in the continuum. 

We find the emergent gauging of translational sym- 
metry realized in the quantum nematic an even better 
example of the usefulness of this insight. Earlier work 
indicated that the relativistic version of this nematic is 
somehow associated with emergent gravity. Resting on 
the "coherent superposition" argument it becomes di- 
rectly transparent what causes the gauging of the crys- 
tal coordinates: the condensed dislocations "shake the 
coordinates coherently" such that infinitesimal Einstein 
translations appear while the Lorentz frame stays fixed. 
This emergent symmetry imposes that the collective exci- 
tations of the quantum nematic have to be in one-to-one 
correspondence with linearized gravity. 

Our message is that we have identified a mechanism 
for the "dynamical generation" of gauge symmetry which 
is very simple but also intriguing viewed from a general 
physics perspective: the quantum mechanics principle of 
states in coherent superposition being "equally true at 
the same time" translates directly to the principle that 
the global symmetry that is broken in the ordered state 
is turned into a gauge symmetry on the disordered side 
just by the quantum un-determinedness of the topolog- 
ical excitations in the dual condensate. This raises the 
interesting question: is quantum coherence required for 
the emergence of local symmetry, or can it also occur in 
classical systems? 

This question relates directly to the spectacular dis- 
covery of "Dirac monopoles" in spin ice^S. Castelnovo, 
Moessner and Sondhi 47 realized that the manifold of 



ground states ( "frustration volume" ) of this classical geo- 
metrically frustrated spin problem is coincident with the 
gauge volume of a compact U(l) gauge theory, with the 
ramification that it carries Dirac monopoles as topolog- 
ical excitations. All along it has been subject of debate 
to what extent these monopoles can be viewed as lit- 
eral Dirac monopoles in the special "vacuum" realized in 
the spin ice, or rather half-bred cartoon versions of the 
real thing. With our recipe at hand it is obvious how 
to make them completely real: imagine the classical spin 
ice to fill up Euclidean spacetime, and after Wick rota- 
tion our "coherent superposition principle" would have 
turned the frustration volume of the classical problem 
into a genuine gauge volume since by quantum superpo- 
sition all degenerate states would be "equally true at the 
same time" . 

The ambiguity associated with the classical spin ice 
monopoles is rooted in the role of time. In principle, by 
doing time-resolved measurements one can observe every 
particular state in the frustration volume and this ren- 
ders these states to be not gauge equivalent. However, all 
experiments which have revealed the monopoles involved 
large, macroscopic time scales. One can pose the ques- 
tion whether it is actually possible under these conditions 
to define observables that can discriminate between the 
"fake" monopoles of spin ice and the monopoles of Dirac. 
Perhaps the answer is pragmatic: as long as ergodicity 
is in charge, one can rely on the ensemble average in- 
stead of the time average, and as long as the time scale 
of the experiment is long enough such that one is in the 
ergodic regime, the frustration volume will "disappear" 
in the ensemble average. For all practical purposes one is 
then dealing with a genuinely emergent gauge symmetry 
which tells us that in every regard the spin ice monopole 
is indistinguishable from the Dirac monopole. 
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